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potential-energy  expressions  have  recently  been  proposed 
for  carbon,  11  no  potential-energy  functions  are 
currently  available  that  can  be  used  in  molecular- 
dynamics  simulations  to  model  the  CVD  of  diamond 
films  from  hydrocarbons. 

As  a  first  step  toward  using  molecular  dynamics  to 
model  the  CVD  of  diamond  films,  we  have  developed  an 
empirical  potential-energy  function  that  captures  many 
of  the  essential  features  of  chemical  bonding  in  hydrocar¬ 
bons.  In  particular,  the  goal  of  this  work  has  been  to  de¬ 
velop  a  simple  expression  that  (1)  can  reproduce  the  in¬ 
tramolecular  energetics  and  bonding  in  solid  diamond 
and  graphite  as  well  as  a  number  of  essential  hydrocar¬ 
bon  molecules;  (2)  yields  realistic  energetics  and  bonding 
for  structures  not  included  in  the  data  base;  (3)  allows  for 
bond  breaking  and  forming  (i.e.,  chemistry);  and  (4)  is  not 
computationally  intensive.  While  a  number  of  classical 
expressions  have  been  developed  for  studying  the  dynam¬ 
ics  of  hydrocarbon  molecules,  none  are  able  to  meet  these 
criteria.  For  example,  valence  force  fields  and  related  ap¬ 
proaches  such  as  the  molecular  mechanics  formalism 
developed  by  Allinger  and  co-workers12  do  not  incorpo¬ 
rate  bond  dissociation,  and  so  cannot  be  used  to  model 
the  chemical  addition  of  molecules  to  diamond  surfaces. 
The  approach  taken  here  is  to  use  a  highly  parametrized 
Tersoff  bond-order  expression11  as  a  fitting  function  for  a 
number  of  small  molecules  as  well  as  graphite  and  dia¬ 
mond  lattices.  Since  a  bond-order  function  is  fit  rather 
than,  for  example,  a  many-body  expansion  of  the  poten¬ 
tial  energy,  the  function  should  be  reasonably  transfer¬ 
able  to  structures  not  included  in  the  fitting  data  set. 

The  organization  of  this  paper  is  as  follows.  In  the 
next  section  the  Abell-Tersoff  bonding  formalism  is  out¬ 
lined,  and  several  problems  associated  with  describing 
carbon  and  hydrocarbons  are  discussed.  Two  parame- 
trizations  of  a  Tersoff-type  expression  are  then  given 
which  have  been  fit  to  the  energetics  of  a  relatively  large 
number  of  hydrocarbon  molecules.  The  first  parametri- 
zation  provides  a  good  description  of  carbon-carbon  bond 
lengths  but  has  stretching  force  constants  that  are  too 
small,  while  the  second  yields  force  constants  that  are  in 
better  agreement  with  experiment  but  produces  bond 
lengths  for  double  and  triple  bonds  that  are  too  long. 
Transferability  of  the  expressions  is  tested  by  calculating 
energetics  for  diamond  surfaces  and  hydrocarbons  not 
used  in  the  fitting  data  base.  The  energetics  of  a  limited 
number  of  atomic  configurations  that  are  of  direct 
relevance  to  the  CVD  of  diamond  films  is  then  given. 

II.  POTENTIAL-ENERGY  EXPRESSION 

In  an  attempt  to  explain  a  universality  relation  ob¬ 
served  in  binding-energy  curves,14  Abell  derived  a  general 
expression  for  binding  energy  that  is  a  sum  of  near¬ 
neighbor  pair  interactions  that  are  moderated  by  the  lo¬ 
cal  atomic  environment.15  Tersoff  subsequently  intro¬ 
duced  an  analytic  potential-energy  function  based  on  the 
Abell  expression  that  realistically  describes  bonding  in 
silicon  for  a  large  number  of  solid-state  structures.11 
Similar  expressions  have  since  been  used  for  a  range  of 
applications,  including  other  IV-IV  (Refs.  28,  29,  36,  and 


37)  and  III-V  compounds,11  as  a  fitting  function  for  the 
reactions  H+H2  and  O+O2,29  and  as  a  model  for 
describing  reactive  collisions  in  molecular  solids.19 

The  binding  energy  in  the  Abell-Tersoff  formalism  is 
written  as  a  sum  over  atomic  sites  t, 

=  7  «) 

/ 

where  each  contribution  Ej  is  written  as 

£,=  2I V*ir,j)-B,jVAlru)]  .  (2) 

jt*n 

In  Eq.  (2),  the  sum  is  over  nearest  neighbors  j  of  atom  /, 
yR(r)  and  VA(r)  are  pair-additive  repulsive  and  attrac¬ 
tive  interactions,  respectively,  and  By  represents  a 
many-body  coupling  between  the  bond  from  atom  1  to 
atom  j  and  the  local  environment  of  atom  i.  As  discussed 
by  Abell  and  Tersoff,11, 15  if  Morse-type  functions  are  used 
for  the  repulsive  and  attractive  pair  terms,  then  By  can 
be  considered  a  normalized  bond  order  because  the  Paul¬ 
ing  relationship  between  bond  order  and  bond  length  is 
realized.40 

Abell  has  suggested  that  to  a  first  approximation  By 
can  be  given  as  a  function  of  local  coordination  Z  in  the 
form 

By*Z'h ,  (3) 

where  8  may  depend  on  the  particular  system.  -This 
analysis  yields  a  relationship  between  bond  length,  bind¬ 
ing  energy,  and  coordination  that  is  realized  for  a  number 
of  systems,  including  carbon-carbon  bonds  in  diamond, 
graphite,  and  in  hydrocarbon  molecules.27 

The  procedure  used  by  Tersoff  to  develop  classical  po¬ 
tentials  for  silicon,  carbon,  and  germanium  has  been  to  fit 
the  pair  terms  and  an  analytic  expression  for  By  to  a 
number  of  properties  of  the  diatomic  and  solid-state 
structures  (e.g.,  bond  energies  and  lengths,  bulk  moduli, 
vacancy  formation  energies,  etc.).  His  expression 
developed  in  this  way  appears  to  be  relatively  transferable 
to  other  solid-state  structures  not  used  in  the  fitting  pro¬ 
cedure  such  as  surface  reconstructions  on  silicon11  and 
interstitial  defects  in  carbon.28  As  discussed  below,  how¬ 
ever,  further  analysis  shows  that  this  expression  is  unable 
to  reproduce  a  number  of  properties  of  carbon  such  as  a 
proper  description  of  radicals  and  conjugated  versus  non- 
conjugated  double  bonds. 

A.  7-bonding  in  carbon 

Since  the  Abell-Tersoff  expression  uses  Pauling  bond- 
order  ideas,  it  can  realistically  describe  carbon-carbon 
single,  double,  and  triple  bond  lengths  and  energies  in  hy¬ 
drocarbons  and  in  solid  graphite  and  diamond.  In  inter¬ 
mediate  bonding  situations,  however,  the  assumption  of 
near-neighbor  interactions  combined  with  the  sum  over 
atomic  sites  [Eq.  (1)]  results  in  nonphysical  behavior.  For 
example,  if  a  carbon  atom  with  three  nearest  neighbors  is 
bonded  to  a  carbon  atom  with  four  neighbors,  Eqs. 
(l)-(3)  interpolate  the  bond  so  that  it  is  intermediate  be¬ 
tween  a  single  and  double  bond.  The  formation  of  a  dou¬ 
ble  bond,  however,  results  from  the  overlap  of  unbonded 
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2 p  orbitals.  Since  the  atom  with  a  coordination  of  4  does 
not  have  a  free  orbital,  v  overlap  cannot  occur  and  the 
bond  is  better  described  as  a  single  bond  plus  a  radical  or¬ 
bital.  This  overbinding  of  radicals  results  in  a  nonphysi¬ 
cal  description  of  bonding  for  carbon  in  a  number  of 
common  situations.  For  example,  the  formation  of  a  va¬ 
cancy  in  diamond  results  in  four  radicals,  and  so  this 
overbinding  of  radicals  makes  fitting  the  vacancy  forma¬ 
tion  energy  in  diamond  while  maintaining  a  fit  to  graph¬ 
ite  difficult  if  not  impossible.21 

Similar  nonphysical  behavior  results  when  conjugated 
and  nonconjugated  double  bonds  are  examined.  For  ex¬ 
ample,  in  graphite  each  atom  has  a  local  coordination  of 
3,  and  Tersoff-type  potentials  for  carbon  have  been  fit  to 
yield  a  bond  strength  appropriate  for  graphite  for  this 
atomic  environment.28,2’  A  simple  analysis  of  kekule’ 
structures  for  graphite  shows  that  due  to  the  conjugation 
each  bond  has  approximately  one-third  double-bond  and 
two-thirds  single  bond  character.  In  the  molecule 
(CHj)2C=C(CHj)2  the  two  carbon  atoms  connected  by 
the  central  bond  have  the  fame  local  environment  as  in 
graphite,  but  because  the  bond  is  not  conjugated  it  has  an 
almost  entirely  dot  1  le-bond  character.  Hence  the  poten¬ 
tial  cannot  describe  both  situations  unless  nonlocal  effects 
are  included. 

One  way  to  correct  both  of  these  problems  while  main¬ 
taining  the  fit  to  diamond  and  graphite  is  to  rewrite  Eqs. 
(1)— (3)  as  a  sum  over  bonds  in  the  form 

s*-2  2  IWv (4) 

/  JOi) 

where 

Bij={Blj+BJi)/2  .  (5) 

The  overbinding  of  radicals  can  now  be  fixed  by  adding 
corrections  to  Eq.  (5)  for  bonds  between  pairs  of  atoms 
that  have  different  coordinations.  As  described  below, 
nonlocal  effects  can  also  be  added  to  a  first  approxima¬ 
tion  to  account  for  conjugated  versus  nonconjugated 
bonding. 


B.  Empirical  potential-energy  expression  for  hydrocarbons 

Following  the  discussion  given  above,  the  binding  ener¬ 
gy  for  the  hydrocarbon  potential  is  given  as  a  sum  over 
bonds  as 

^“2  2  1^*)-^^)].  (6> 

I  jl>il 

where  the  repulsive  and  attractive  pair  terms  are  given  by 

)D}j*/(Sjj  - 1  (7) 

and 

)D}f'S,jAS,j  - 1  k ,  (8) 

respectively.  The  function  /y(r),  which  restricts  the  pair 
potential  to  nearest  neighbors,  is  given  by 

1,  r<Rlj" 

fij(r)=  1+cos 

0,  r>Riij)  . 

(9) 

The  pair  terms  are  of  the  same  form  as  that  used  by  Ter- 
joflr28.33.36  jjUt  rewritten  to  make  the  correspondence  to  a 
Morse  function  more  apparent.  If  SfJ  — 2,  then  the  pair 
terms  reduce  to  the  usual  Morse  potentiaL  Furthermore, 
the  well  depth  Djf  \  equilibrium  distance  R\p  and  P,j  are 
equal  to  the  usual  Morse  parameters  independent  of  the 
value  of  Sy. 

The  empirical  bond-order  function  is  given  by  the  aver¬ 
age  of  terms  associated  with  each  atom  in  a  bond  plus  a 
correction  as  discussed  above: 

B^iBij+Bji )/2+Fij(Nj'\ Nj'\ )  ,  (10) 

where 


irif-Rj")  / 


Riy,<r<R^» 


1+  2  ( Oijk  )fik ( riit  )e a‘il [< r<y  <r* 

k(*i.j) 


The  quantities  NjC)  and  NjH)  are  the  number  of  carbon 
and  hydrogen  atoms,  respectively,  bonded  to  atom  /,  Nj'] 
is  the  total  number  of  neighbors  (NjC)+ NjH))  of  atom  i, 
Nfj0”1  depends  on  whether  a  bond  between  carbon  atoms  i 
and  j  is  part  of  a  conjugated  system,  G ( 6,Jk )  is  a  function 
of  the  angle  between  bonds  i — j  and  / — k,  and  the  func¬ 
tions  H,j  and  Ftj  are  described  below.  The  latter  func¬ 
tion,  which  is  used  for  carbon-carbon  bonds  only,  is  the 
correction  discussed  above.  Forms  similar  to  that  used 
by  Tersoff  to  describe  group-IV  alloys  were  tried,26  but 
because  of  the  large  difference  in  bonding  characteristics 
between  hydrogen  and  carbon  (hydrogen  is  monovalent 
while  carbon  has  a  valency  of  up  to  4),  a  set  of  parame¬ 


ters  could  not  be  found  that  could  adequately  describe 
bond  energies  for  a  large  number  of  hydrocarbon  mole¬ 
cules.  Furthermore,  as  discussed  above,  radicals  and 
nonconjugated  double  ly  nds  are  not  well  described. 

To  make  the  potential  function  continuous,  the  func¬ 
tions  are  used  to  define  bonding  connectivity  in  the 
system,  and  are  therefore  used  to  define  values  for 
Nja,  Nj'\  and  Nff'K  First,  values  for  NjH\  NjC),  and 
]V(,rt  for  each  of  the  carbon  atoms  i  are  given  by 

NjH)=  2  ftj(rjj)  ,  (12) 

ji  m hydrogen) 

*;c,=  2  w-  (13) 

ji  “carbon) 


42 


EMPIRICAL  POTENTIAL  FOR  HYDROCARBONS  FOR  USE  IN  . . . 


9461 


and 

N}''=NjH'+NjC) .  (14) 

Values  of  N}'1  for  neighbors  of  the  two  carbon  atoms  in¬ 
volved  in  a  bond  can  be  used  to  determine  whether  the 
bond  is  part  of  a  conjugated  system.  For  example,  if  any 
neighbors  are  carbon  atoms  that  have  a  coordination  of 
less  than  4  (JV/'1  <4),  the  bond  is  defined  as  being  part  of 
a  conjugated  system.  The  value  of  Ny*  for  a  bond  be¬ 
tween  carbon  atoms  i  and  j  is  given  by 

JV/r=  1+  2  fikirik)F(xik) 

cartoons  Jkt 

+  2  f/r^FUj,)  (15) 

carbons  H*iJ) 


energies.44  For  methane,  the  heat  of  formation  at  0  K  of 
0.693  eV  (Ref.  41)  was  used  to  determine  carbon- 
hydrogen  bond  energies  of  4.393  eV.  Bond  dissociation 
energies  of  4.861  eV  for  methane45,  and  4.944  eV  for 
methyl  radical46  were  then  used  to  determine  bond  ener¬ 
gies  of  4.237  eV  for  CH3  and  3.833  eV  for  CH2. 

Two  different  parametrizations  for  Eqs.  (6)— (18)  for 
carbon  and  hydrocarbons  have  been  determined.  A  simi¬ 
lar  function  (without  the  corrections  discussed  above)  has 
been  fit  to  the  binding  energy  of  the  Cj  diatomic  mole¬ 
cule,  and  the  binding  energies  and  lattice  constants  of 
graphite,  diamond,  simple  cubic,  and  face-centered-cubic 
structures,29  and  so  the  first  parametrization  was  chosen 
to  reduce  to  that  function.  The  function  Gc(0)  is  given 
by 


where 


FUik)= 


1,  xik<2 

( 1  +cos[ idxik—  2)1} /2,  2<xik<2  (16) 

0,  xik  >  3 


and 

x,k=Nr-fiklrik)  ■  (17) 

This  function  yields  a  continuous  value  of  N con*  as  bonds 
break  and  form  and  as  second-neighbor  coordinations 
change.  For  lV“nj=  1,  a  bond  is  not  part  of  a  conjugated 
system  and  the  function  yields  appropriate  values,  and 
for  2  the  bond  is  part  of  a  conjugated  system  and 

parameters  fit  to  conjugated  bonds  are  used.  As  a  final 
step  for  making  the  potential  continuous,  two-  and 
three-dimensional  cubic  splines  are  used  for  the  functions 
H,j  and  FtJ,  respectively,  to  interpolate  between  values  at 
discrete  numbers  of  neighbors. 

The  procedure  used  to  determine  parameters  for  Eqs. 
(6)— (17)  was  to  first  fit  to  systems  consisting  of  only  car¬ 
bon  and  only  hydrogen.  Parameters  were  then  chosen  for 
the  mixed  hydrocarbon  system  that  reproduced  additive 
bond  energies.  Since  the  pair  terms  are  first  fit  to  solid- 
state  carbon  structures,  the  equilibrium  carbon-carbon 
distances  and  stretching  force  constants  for  hydrocarbons 
are  completely  determined  by  fitting  to  bond  energies. 

To  determine  appropriate  energies  for  hydrocarbons 
with  carbon-carbon  bonds,  additive  bond  energies  for  sin¬ 
gle,  double,  conjugated  double  and  triple  carbon-carbon 
bonds,  and  carbon-hydrogen  bonds  were  determined 
from  molecular  atomization  energies.  Heats  of  formation 
at  0  K  for  acetylene,  ethylene,  ethane,  benzene,  and  cy¬ 
clohexane  are  2.356,  0.629,  —0.717,  1.041,  and  —0.868 
eV,  respectively.41  Combining  these  values  with  binding 
energies  of  7.3768  eV  for  graphite42  and  2.375  eV  for 
molecular  hydrogen,45  and  neglecting  corrections  for  zero 
point  energy  in  the  hydrocarbons,  yields  total  atomiza¬ 
tion  energies  of  17.149,  23.626,  29.723,  57.472,  and 
73.634  eV  for  acetylene,  ethylene,  ethane,  benzene,  and 
cyclohexane,  respectively.  Assuming  constant  carbon- 
hydrogen  and  carbon-carbon  single  bond  energies,  energy 
values  of  4.362,  8.424,  6.175,  5.216,  and  3.547  eV  can  be 
derived  for  carbon-hydrogen  and  carbon-carbon  triple, 
double,  conjugated  double  (in  benzene),  and  single  bond 


Gc(0)=flo|l+cS/4-cM4+<l+cos0)2])  «« 

and  the  carbon-carbon  parameters  are  given  in  Table  I. 
The  potential  yields  binding  energies  and  nearest- 
neighbor  bond  lengths  of  7.3768  eV/atom  and  1.42  A,  re¬ 
spectively,  for  graphite,  and  7.3464  eV/atom  and  1.54  A 
for  diamond.  It  also  reproduces  calculated  values  of  7.2 
and  7.6  eV  for  the  vacancy  formation  energies  in  dia¬ 
mond  and  graphite,  respectively.4*  The  predicted 
carbon-carbon  bond  lengths  and  stretching  force  con¬ 
stants  given  by  this  parametrization  (denoted  as  potential 
I)  for  the  bond  energies  derived  above  are  given  in  Table 
II.  The  bond  lengths  agree  very  well  with  the  expenmen- 
tal  values,  the  largest  difference  being  0.02  A  for  a  single 
bond.  The  stretching  force  constants,  however,  differ 
from  experimental  values  by  as  much  as  60%.  The 
carbon-hydrogen  bond  length  given  by  the  potential  is 
1.071  A,  which  is  slightly  short  of  typical  experime.  ‘al 
values  of  1.09  A.40 

Given  in  Table  III  are  a  second  set  of  hydrocarbon  pa¬ 
rameters  for  Eqs.  (6)- (18).  This  potential  was  again  fit  to 
the  binding  energies  and  lattice  constants  of  graphite,  di¬ 
amond,  simple  cubic,  and  face-centered-cubic  structures 
for  pure  carbon,42  and  the  vacancy  formation  energies  for 
diamond  and  graphite.4*  It  yields  binding  energies  and 
nearest-neighbor  bond  lengths  of  7.3756  eV/atom  and 
1.45  A,  respectively,  for  graphite,  and  7.3232  eV/atom 
and  1.54  A  for  diamond.  The  predicted  carbon-carbon 
bond  lengths  and  force  constants  for  hydrocarbons  given 
by  this  parametrization  (denoted  by  potential  II)  are 
given  in  Table  II.  The  stretching  force  constants  are 
closer  to  the  experimental  values  than  potential  I,  with 
the  differences  being  reduced  to  11%,  8%,  and  26%  for 
single,  double,  and  triple  bonds,  respectively.  The  bond 
lengths  for  multiple  bonds,  however,  are  not  as  well  de¬ 
scribed  as  with  potential  I,  with  the  double  and  triple 
bonds  being  larger  than  the  experimental  values  by  3.7% 
and  7.5%,  respectively.  A  parametrization  for  the  func¬ 
tions  used  for  the  pair  terms  could  not  be  found  that  fit 
carbon-carbon  stretching  force  constants  and  bond 
lengths  simultaneously,  although  other  forms  may  repro¬ 
duce  both  sets  of  properties. 

For  hydrogen,  the  pair  terms  were  fit  to  properties  of 
the  gas-phase  diatomic.45  The  parameter  SH  was  set 
equal  to  8C,  and  the  remaining  parameters  were  fit  to 
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TABLE  I.  Parameters  for  Eqs.  (4)-(18)  for  potential  I.  All  parameters  not  given  are  equal  to  zero,  F(i,j,k)=* FI j,i,k), 
U,j,k  >2)—F(i,j,2),  and  dFU,j,k)/di=dF(j,i,k)/di.  The  partial  derivatives  are  used  in  the  multidimensional  cubic  splines. 


Carbon 


Hydrogen 


rameter 

Value 

Fit  to 

Diet 

Acc 

1.315  A 

Lattice  constants  of  diamond  and  graphite 

UCC 

6.325  eV 

Reference  29 

Pcc 

1.5  A-' 

Reference  29 

Sec 

1.29 

Reference  29 

8Cc 

0.80469 

Reference  29 

accc 

0.0 

pill 

Acc 

1.7  A 

Reference  29 

P’121 

«cc 

2.0  A 

Reference  29 

0.011304 

Reference  29 

Co 

19J 

Reference  29 

dl 

2.5* 

Reference  29 

F(2,3,i),F(2,3,2) 

-0.0465 

for  diamond 

F(  1,2,2) 

-0.0355 

En c  for  graphite 

p  Irl 
*HH 

0.74144  A 

Gas-phase  diatomic 

r\ie) 

4.7509  eV 

Gas-phase  diatomic 

Phh 

1.9436  A'1 

Gas-phase  diatomic 

$hh 

2.3432 

Barrier  for  reaction  (19) 

Shh 

0.80469 

Set  equal  to  carbon  value 

aHHH 

3.0 

Remove  spurious  wells  from  (19) 

CHH 

4.0 

Barrier  for  reaction  (19) 

P(ll 
**  HH 

1.1  A' 

Near-neighbor  interactions 

p(2l 
**  HH 

1.7  A 

Near-neighbor  interactions 

p  If! 

ach 

1.1199  A 

Gas-phase  diatomic 

p\le» 

^CH 

3.6422  eV 

Gas-phase  diatomic 

Pdt 

1.9583  A*1 

Gas-phase  diatomic 

1.7386 

i^HH^cc),/J 

nil) 

**  CH 

1.3  A 

Reaction  (19) 

p(2» 

aCH 

1.8  A 

Reaction  (19) 

aHHOaCHH»aHCH>aHCC 

3.0  A-1 

Hj  value 

Hcc(l,l) 

-0.0175 

CC  bond  energy  in  benzene 

//cc(2,0) 

-0.0070 

CC  double  bond  in  ethylene 

Wcc(3,0) 

0.0119 

CC  single  bond  in  ethane 

//ccU.2) 

0.0115 

CC  single  bond  in  isobutane 

«cc(2.1) 

0.0118 

CC  single  bond  in  cyclohexane 

HchU.O) 

-0.0760 

Atomization  energy  of  CH2 

«cm(2,0) 

-0.2163 

Atomization  energy  of  CHj 

«ch(3,0) 

-0.3375 

Atomization  energy  of  methane 

■^CHiO>l) 

-0.1792 

CH  bond  energy  in  acetylene 

tfCH<0,2) 

-0.2407 

CH  bond  energy  in  benzene 

«ch(1.1) 

-0.2477 

CH  bond  energy  in  ethylene 
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TABLE  I.  <  Continued). 


Parameter 

Value 

Fit  to 

-0.3320 

CH  bond  energy  in  ethane 

-0.3323 

Tertiary-HC  bond  energy  in  isobutane 

Hc *0.2) 

-0.3321 

CH  bond  energy  in  cyclohexane 

3Hc«(l.l) 

dC 

-0.12805 

Centered  difference 

dHcn.o) 

ac 

-0.07655 

Centered  difference 

a  i/c««U) 

3  H 

-0.13075 

Centered  difference 

-0.0764 

Centered  difference 

3  H 

- 

FU.1,1) 

0.1511 

CC  triple  bond  in  acetylene 

R2,2,l> 

0.075 

Average  energy  of  bonds  in  (CHj)2C=C(CHj) 
and  (CHj)HC=CH(CHj)  equal  double  bond 

FU.2,1) 

0.0126 

Atomization  energy  of  HC=CHj 

FlU.D^lU^) 

-0.1130 

Single  bond  in  H3C — CH 

FTO,3,l),m3,2) 

-0.1220 

Single  bond  in  HjC — C 

FI0A2) 

-0.0445 

Conjugated  double  bond  in  C=CH(CH2) 

Ft  0,2,1) 

0.0320 

Double  bond  in  C— CH2 

FlO.l.l) 

0.1100 

Atomization  energy  of  C2H 

fU.1,2) 

0.0074 

Atomization  energy  of  CHjCCH 

3F(3, 1,1) 

-0.1160 

Centered  difference 

3/ 

3F<3,2,1) 

3i 

-0.13205 

Centered  difference 

3F(3,1,2) 

3/ 

-0.0610 

Centered  difference 

3F(2,3,2) 

3i 

0.02225 

Centered  difference 

3F(2,4,2) 

3  i 

-0.037  75 

Centered  difference 

3F(3,4,2) 

3 1 

0.0565 

Centered  difference 

3F(3.4.I) 

3/ 

0.0565 

Centered  difference 

3F<3,2,2) 

-0.1065 

Centered  difference 

3/ 

reproduce  the  barrier  of  0.425  eV  for  the  ground-state 
linear  exchange  reaction 

H+Hj-*H2+H  (19) 

at  the  hydrogen-hydrogen  distance  of  0.93  A.  as  calculat¬ 
ed  by  Liu  and  Seigbahn49  and  to  remove  spurious  wells 
from  this  surface.  For  simplicity,  angular  dependence  of 
the  potential  for  hydrogen  centers  via  the  function  GH(9) 
was  replaced  by  the  (constant)  value  of  GH  used  to  fit  re¬ 
action  (19).  The  resulting  parameters  corresponding  to 


hydrocarbon  potential  I  are  given  in  Table  I,  and  those 
for  potential  II  are  given  in  Table  III. 

Some  important  points  about  the  function,  fitting  pro¬ 
cedure,  and  resulting  structures  need  to  be  addressed. 
First,  the  angular  function  G(0)  associated  with  the  car¬ 
bon  centers  favors  180*  bond  angles  and  hence  open 
structures.  This  is  physically  motivated  by  valence  shell 
electron  pair  repulsion  (VSEPR)  theory,50  which  assumes 
that  repulsions  between  pairs  of  valence  electrons  tend  to 
maximize  bond  angles.  This  assumption,  combined  with 
the  tendency  to  maximize  the  binding  energy  (i.e.,  num- 
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TABLE  II.  Carbon-carbon  bond  lengths  R,  (in  A)  and  force  constants  Fc  (in  10*  dyn/cm)  given  by 
the  two  potentials  based  on  bond  energies  of  3.547,  6.175,  and  8.424  eV  for  single,  double,  and  triple 
bonds,  respectively,  as  derived  in  the  text. _ 


Single  bonds 

Fe  R, 

Double  bonds 

Ft  Rr 

Triple  bonds 

Fc  R, 

Potential  I 

2.6 

1.56 

4.5 

1.33 

6.1 

1.20 

Potential  II 

5.0 

1.55 

8.7 

1.38 

11.9 

1.29 

Experiment 

4.5* 

1.54k 

9.5* 

1.33" 

16.0* 

l^O" 

•From  Ref.  40. 
bFrom  Ref.  48. 


TABLE  III.  Parameters  for  Eqs.  (4)— ( 18)  for  potential  II.  All  parameters  not  given  are  equal  to  zero,  F{i,j,k)~ F(j,i,k), 
F(.i,j,k  >  2)=F(/,y,2),  and  dF(i,j,k)/di=dF{j,i,k) /3i.  The  partial  derivatives  are  used  in  the  multidimensional  cubic  splines. _ 


Parameter 

Value 

Fit  to 

Carbon  R'cc 

1.39  A 

Lattice  constants  of  diamond  and  graphite 

r\ie) 

UCC 

6.0  eV 

Fit  to  data  in  Ref.  42 

Pcc 

LI  A-1 

Fit  to  data  in  Ref.  42 

•Sec 

1.22 

Fit  to  data  in  Ref.  42 

Soc 

0.5 

Fit  to  data  in  Ref.  42 

«ccc 

0.0 

nil) 

KCC 

1.7  A 

Fit  to  data  in  Ref.  42 

n  (2) 

^CC 

2.0  A 

Fit  to  data  in  Ref.  42 

ao 

0.00020813 

Fit  to  data  in  Ref.  42 

Co 

330* 

Fit  to  data  in  Ref.  42 

-  dl 

3.5J 

Fit  to  data  in  Ref.  42 

5123,11,503.2) 

-0.0363 

E„c  for  diamond 

511 A2). 

-0.0243 

for  graphite 

Hydrogen  R'hh 

0.74144  A 

Gas-phase  diatomic 

4.7509  eV 

Gas-phase  diatomic 

0HH 

1.9436  A-' 

Gas-phase  diatomic 

•Shh 

L3432 

Barrier  for  reaction  (19) 

x 

X 

0.5 

Set  equal  to  carbon  value 

aHHH 

4.0 

Remove  spurious  wells  from  (19) 

12.33 

Barrier  for  reaction  (19) 

p«l> 

**  HH 

1.1  A 

Near-neighbor  interactions 

p<2> 

^HH 

1.7  A 

Near-neighbor  interactions 

Hydrocarbons  Rq» 

1.1199  A 

Gas-phase  diatomic 

UCH 

3.6422  eV 

Gas-phase  diatomic 

Pen 

1.9583  A-' 

Gas-phase  diatomic 

$CH 

1.69077 

(ShhScc),/5 

pill 

1.3  A 

Reaction  (19' 

p  12) 

1.8  A 

Reaction  (19) 

H  HC»  Ooih,  aHCH » aHCC 

4.0  A"' 

H,  value 

EMPIRICAL  POTENTIAL  FOR  HYDROCARBONS  FOR  USE  IN  . . . 
TABLE  III.  ( Continued ). 


Parameter 

Value 

Fit  to 

Hcc(  1,1) 

-0.0226 

CC  bond  energy  in  benzene 

Wcc<2.0> 

-0.0061 

CC  double  bond  in  ethylene 

Hcc(  3,0) 

0.0173 

CC  single  bond  in  ethane 

Hcc(\,2) 

0.0149 

CC  single  bond  in  isobutane 

Hcc(  2,1) 

0.0160 

CC  single  bond  in  cyclohexane 

Hch(1,0) 

-0.0984 

Atomization  energy  of  CHi 

Hch(2,0) 

-0.2878 

Atomization  energy  of  CHj 

HchU,0) 

-0.4507 

Atomization  energy  of  methane 

Hc  „<0,1) 

-0.2479 

CH  bond  energy  in  acetylene 

HoitO.2) 

-0.3221 

CH  bond  energy  in  benzene 

Hch(U) 

-0.3344 

CH  bond  energy  in  ethylene 

Hch(2.U 

-0.4438 

CH  bond  energy  in  ethane 

Hch(0,3) 

-0.4460 

Tertiary-HC  bond  energy  in  isobutane 

HCh(1.2» 

-0.4449 

CH  bond  energy  in  cyclohexane 

a^cHiU) 

3  C 

-0.17325 

Centered  difference 

3i/cH<2,0) 

3  C 

-0.09905 

Centered  difference 

3//ch(0,2) 

3tf 

-0.17615 

Centered  difference 

3//ch(1.D 

3  H 

-0.09795 

Centered  difference 

HI, 1,1) 

0.1264 

CC  triple  bond  in  acetylene 

H  2,2,1) 

0.0605 

Average  energy  of  bonds  in  (CHj)jC=C(CHj) 

and  (CHj)HC=CH(CHj)  equal  double  bond 

/IU.II 

0.0120 

Atomization  energy  of  HC=CHj 

HI, 3, 1)^11, 3^) 

-0.0903 

Single  bond  in  HjC — CH 

HO,3,l)H(0,3r2) 

-0.0904 

Single  bond  in  HjC — C 

HO, 2, 2) 

-0.0269 

Conjugated  double  bond  in  C=CH(CH:) 

H0A1) 

0.0427 

Double  bond  in  C=CH2 

HO, 1,1) 

0.0996 

Atomization  energy  of  C2H 

HI, 1,2) 

0.0108 

Atomization  energy  of  CH2CCH 

3F(3, 1,1) 

3/ 

-0.0950 

Centered  difference 

6F(  3,2, 1 ) 

3/ 

-0.108  35 

Centered  difference 

3F(3, 1,2) 

3/ 

-0.0452 

Centered  difference 

3F12,3,2) 

3/ 

3F(2,4,2) 

3/ 

3F(3,4,2) 

3/ 

3F(  3,4,1) 
3j 

3F(3,2,2) 

li 


0.01345 
-0.02705 
0.04515 
0.045  15 


Centered  difference 
Centered  difference 
Centered  difference 
Centered  difference 


-0.087  60 


Centered  difference 
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ber  of  bonds  times  the  bond  strength),  leads  to,  for  exam¬ 
ple,  a  tetrahedral  structure  for  methane  and  linear  struc¬ 
ture  for  acetylene.  Because  lone  electron  pairs  are  not  ex¬ 
plicitly  treated  in  this  empirical  potential,  however,  some 
structures  are  linear  that  experimentally  are  bent,  such  as 
methylene.51  Second,  although  the  procedure  for  defining 
Nfj"*  incorporates  conjugation  in  an  approximate  way,  it 
will  not  give  the  subtle  differences  in  bond  orders  for  con¬ 
jugated  ring  systems  predicted  by  even  simple  molecular- 
orbital  methods.  This  procedure,  however,  is  able  to  in¬ 
clude  to  a  first  approximation  nonlocal  effects  without 
the  need  for  diagonalizing  a  matrix  or  going  beyond 
nearest-neighbor  interactions,  hence  drastically  reducing 
computational  requirements  for  large  systems.  Third,  the 
carbon-carbon  single  bond  energies  derived  for  hydrocar¬ 
bons  are  smaller  than  those  in  diamond.  This  difference 
also  occurs  in  more-elaborate  least-squares-fitting  pro¬ 
cedures  for  determining  carbon-carbon  bond  energies, 
and  can  be  attributed  to  both  the  weakening  of  the  bonds 
due  to  hydrogen  and  steric  repulsions.52  Fourth,  barriers 
for  rotation  about  carbon-carbon  bonds  (especially  dou¬ 
ble  bonds)  and  nonbonded  interactions  such  as  van  der 
Waals  forces  have  not  been  included.  With  the  potential 
written  as  a  sum  over  pairs,  the  former  effects  can  be 
modeled,  and  methods  for  incorporating  both  types  of  in¬ 
teractions  within  this  formalism  are  being  developed.  Fi¬ 
nally,  the  price  one  pays  for  using  a  simple  classical  ex¬ 
pression  is  that  properties  such  as  resonance  effects55  (i.e., 
4n  +2  Huckel  stability)  and  Woodward-Hoffman.  rules55 
may  not  necessarily  be  obeyed.  Care  should  therefore  be 
used  when  interpreting  reaction  mechanisms  suggested 
by  '■  npirical  potentials,  and  where  possible  they  should 
be  reexamined  using  other  less-empirical  techniques. 


IIL  PREDICTED  ENERGETICS  AND  BONDING 

To  test  the  transferability  of  the  potential  given  in  the 
preceding  section,  the  energy  and  structure  of  a  variety  of 
small  hydrocarbon  molecules  and  the  diamond  {111}  and 
(001  j  surfaces  have  been  examined.  Shown  in  Table  IV 
are  atomization  energies  given  by  the  empirical  potential 
and  corresponding  experimental  values.  The  experimen¬ 
tal  values  were  derived  from  heats  of  formation  without 
correcting  for  zero-point  energy  or  other  finite  tempera¬ 
ture  effects.  The  overall  agreement  is  good,  with  energies 
being  reproduced  to  within  1%  or  better  for  81%  of  the 
molecules  tested.  While  this  level  of  agreement  may  be 
somewhat  fortuitous  given  the  approximations  made  in 
estimating  the  experimental  values,  it  does  demonstrate 
that  the  potential  can  describe  chemical  bonding  for 
a  wide  range  of  hydrocarbon  molecules.  The  worst  fits 
are  for  the  molecules  CHjCH=C=CH2  and 
HjCssCsCHj  for  both  potentials,  and  for  cyclopro- 
pene  using  potential  II.  For  the  former  two  molecules 
the  potential  function  treats  the  two  double  bonds  as  part 
of  a  conjugated  system,  and  so  the  bond  energies  are  re¬ 
duced.  This  is  an  example  of  where  our  analytic  approxi¬ 
mation  for  defining  conjugation  breaks  down. 

In  Table  V  the  energetics  predicted  by  the  potential  for 
the  relaxed  bulk-terminated,  Pandey  ir-bonded  chain,54 


and  Chadi  ir-bonded  molecule55  reconstructions  on  the 
diamond  (111)  surface  are  given  along  with  values  calcu¬ 
lated  by  Vanderbilt  and  Louie  using  local-density- 
functional  (LDF)  methods.56  The  values  given  by  the 
empirical  potentials  were  calculated  using  a  slab  ten  lay¬ 
ers  thick  with  16  atoms  per  layer  and  periodic  boundary 
conditions  perpendicular  to  the  surface.  The  bottom  two 
layers  were  held  rigid  and  the  top  eight  layers  were  al¬ 
lowed  to  completely  relax  under  the  influence  of  each  of 
the  potentials  to  the  minimum-energy  configurations 
starting  from  the  positions  reported  by  Chadi  for  the 
various  reconstructions.55  The  potentials  predict  that  the 
(undimerized)  ir-bonded  chain  reconstruction  is  energeti¬ 
cally  preferred,  in  agreement  with  the  LDF  result56  and 
with  a  number  of  experimental  studies.57,51  If  radical  en¬ 
ergetics  is  not  correctly  described,  the  Tersoff-type  poten¬ 
tials  given  in  Refs.  28  and  29  incorrectly  predict  the  re¬ 
laxed  bulk-terminated  surface  to  be  the  most  stable.59 

Experimental  studies  have  shown  that  upon  exposure 
to  atomic  hydrogen  the  ir-bonded  chain  reconstruction 
converts  to  a  hydrogen-terminated  surface  with  the  car¬ 
bon  atoms  reverting  to  positions  that  are  characteristic  of 
the  bulk.57,51,60  To  test  the  potentials,  a  monolayer  of  hy¬ 
drogen  atoms  was  placed  above  the  surface  carbon  atoms 
on  the  relaxed  bulk-terminated  surface  and  the  ir-bonded 
chain  reconstruction  and  the  systems  were  relaxed.  Po¬ 
tential  I  predicts  energies  of  —1.00  eV/(surface  atom) 
and  —0.61  eV/(surface  atom)  for  the  hydrogen-covered 
bulk-terminated  surface  and  hydrogen-covered  ir-bonded 
chain  reconstruction,  respectively,  relative  to  the  clean 
ir-bonded  chain  reconstruction  plus  gas-phase  H2  mole¬ 
cules.  Potential  II  yields  corresponding  energies  of 

—  1.13  eV/(surface  atom)  for  adsorption  on  the  bulk- 
terminated  surface  and  —0.63  eV/(surface  atom)  for  ad¬ 
sorption  on  the  ir-bonded  chain  reconstruction.  Hence 
the  potentials  correctly  predict  that  hydrogen  adsorption 
is  energetically  stable,  and  that  adsorption  on  a  bulk- 
terminated  surface  is  favored  over  the  reconstructed  sur¬ 
face.  Side  views  of  the  equilibrium  structures  for  the  two 
hydrogen-ter.ninated  surfaces  are  shown  in  Fig.  1. 

The  structure  and  energetics  of  the  dimer  reconstruct¬ 
ed  (001 )  surface  of  diamond  are  also  of  interest,  and  de¬ 
tailed  experimental  studies  have  just  recently  been  com¬ 
pleted.61  Potential  I  predicts  a  stabilization  energy  of 

—  5.49  eV/dimer  relative  to  the  bulk-terniinated  surface, 
and  a  surface  dimer  bond  length  of  1.38  A.  Potential  II 
predicts  a  slightly  less  stable  reconstruction  with  an  ener¬ 
gy  of  —5.20  eV/dimer  relative  to  the  bulk-terminated 
surface,  and  a  surface  dimer  bond  length  of  1.43  A. 
These  values  were  obtained  using  a  slab  eight  layers  thick 
with  32  atoms  per  layer.  Periodic  boundaries  were  used 
for  the  two  directions  perpendicular  to  the  surface,  the 
bottom  two  layers  were  held  rigid,  and  the  remainder  of 
the  system  was  relaxed  to  the  minimum-energy  structure 
starting  from  a  surface  reconstruction  consisting  of  a 
2X1  arrangement  of  surface  dimers.  The  dimer  length  is 
stretched  from  the  double  bond  distance  in  ethylene  for 
both  potentials,  and  the  bond  energy  is  correspondingly 
weakened.  This  is  due  to  the  strain  induced  by  the  lat¬ 
tice.  The  monohydride  phase,  where  a  hydrogen  atom  is 
bonded  to  each  atom  in  a  dimer  pair,  was  also  examined 
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TABLE  IV.  Atomization  energies  for  various  hydrocarbon  molecules.  Experimental  values  were  derived  from  heats  of  formation 
using  energies  of  7.3768  cV  for  carbon  and  2.375  eV  for  hydrogen. 


Molecule 

Potential  I 
(eV) 

Potential  II 
(eV) 

Experimental  value 
(eV) 

Alkanes  methane 

17.6 

17.6 

17.6* 

ethane 

29.7 

29.7 

29.7* 

propane 

42.0 

42.0 

42.0* 

n-butane 

54.3 

54.3 

54.3* 

/-butane 

54.3 

54.3 

54.4* 

n-pentane 

66.5 

66.5 

66.6* 

isopentane 

66.5 

66.5 

66.6* 

neopentane 

66.8 

£6.8 

66.7* 

cyclopropane 

35.5 

35.0 

35.8* 

cyclobutane 

48.7 

48.5 

48.2* 

cyclopentane 

61.4 

61.3 

61.4* 

cyclohexane 

73.6 

73.6 

73.6* 

Alkenes  ethylene 

23.6 

23.6 

23.6* 

propene 

36.2 

36.2 

36.0* 

1 -butene 

48.5 

48.5 

48.5b 

cis-butene 

48.8 

48.9 

48.6b 

isobutene 

48.4 

48.4 

48.7b 

(CH3)2C-=C(CH3)j 

73.2 

73.3 

73.4b 

cyclopropene 

28.2 

27.3 

28. 8b 

cyclobutene 

42.4 

42.0 

42.4b 

cyclopentene 

55.7 

55.7 

55.6b 

1,4-pentadiene 

55.0 

55.0 

54. 8b 

CH  j — CHCH  =  CHj 

41.8 

41.9 

42.6b 

CHjCH=C=CHj 

40.4 

40.5 

42.  lb 

HjC=C=CHi 

27.8 

27.9 

29.6b 

Alkynes  acetylene 

17.1 

17  1 

17.1* 

propyne 

29.4 

29.4 

29.7*“ 

1-butyne 

41.7 

41.7 

42.0b 

2-butyne 

41.7 

41.7 

42.2b 

Aromatics  benzene 

57.5 

57.5 

57.5* 

toulene 

69.6 

69.6 

70.  lb 

1 ,4-dimethylbenzene 

81,8 

81.8 

82. 6b 

ethylbenzene 

81.9 

81.9 

82.5b 

ethenylbenzene 

76.2 

76.2 

76.5b 

ethynylbenzene 

69.8 

69.8 

69.9b 

naphthalene 

91.4 

91.4 

91.2b 

Radicals  CHj 

7.8 

7.8 

oo 

A 

CHj 

-  ,  12.7 

12.7 

12.7" 

HAHj 

25.7 

25.7 

25. 5b 

HjCjH 

18.9 

18.9 

18.9* 

C2H 

12.2 

12.2 

12.2f 

CHjCCH 

24.5 

24.5 

25.8b 

n-CjH, 

37.9 

38.0 

37.8' 

i-CjHj 

38.3 

38.3 

38.0* 

t-C,H, 

50.5 

50.5 

50.5' 

phenyl 

52.7 

52.7 

52 .7* 

•From  heat  of  formation  at  0  K  Fro.:,  Tef.  41). 

"From  heat  of  formation  at  30C  K.  (from  H.  M.  Rosenstock,  K.  Draxl,  B.  W.  Steiner,  and  J.  J.  Heron,  J.  Chem.  Ref.  Data  6.  Suppl.  1, 
1-774(1977)]. 

'Calculated  from  P.ef.  46. 

''Calculated  from  Ref.  45. 

'From  heat  of  formation  at  300  K  (from  Ref.  52). 

'From  heat  of  formation  at  0  K  [from  H.  M.  Rosenstock,  K.  Draxl,  B.  W.  Steiner,  and  J.  J.  Heron,  J.  Chem.  Ref.  Data  6,  Suppl.  1, I- 
774(1977)]. 
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TABLE  V.  Energetics  (in  eV/surface  atom)  for  various  structures  of  the  f  1 1 1 1  surface  of  diamond. 


Ideal  1X1 

Relaxed  1 X  1 

tr-bonded 

chain 

ir-bonded 

molecule 

Potential  I 

0.0 

-0.24 

-1.10 

-0.32 

Potential  II 

0.0 

-0.23 

-1.03 

-0.25 

LDP 

0.0 

-0.37 

-0.68 

0.28 

'From  Ref.  56. 

using  the  potentials.  The  predicted  carbon-carbon  sur¬ 
face  dimer  bond  length  is  1.63  A  for  potential  I  and  1.60 
A  for  potential  II,  both  of  which  are  stretched  from  the 
gas-phase  single  bond  lengths.  Potentials  I  and  II  predict 
changes  in  energy  for  hydrogen  adsorption  of  —2.72 
eV/dimer  and  —2.96  eV/dimer,  respectively,  rejative  to 
the  clean  dimer  reconstructed  surface  and  gas-phase  H2. 
Hydrogen  adsorption  is  therefore  predicted  to  be  energet¬ 
ically  stable,  in  agreement  with  experiment.61 

The  abstraction  of  hydrogen  from  hydrogen- 
terminated  diamond  surfaces  has  been  proposed  as  being 
a  rate-determining  step  in  the  CVD  of  diamond  films,1  “3 
so  the  potential-energy  surface  for  the  hydrogen  abstrac- 


FIG.  t.  Side  views  of  hydrogen-terminated  diamond  |111| 
surfaces.  The  top  picture  is  the  rr-bonded  chain  reconstruction 
and  the  bottom  picture  is  the  bulk-terminated  surface. 


tion  reaction 

CH4+H— HjC-H-H— CHj+H2  (20) 

was  examined  using  the  empirical  potentials.  This  reac¬ 
tion  was  chosen  because  the  potential-energy  surface  has 
been  well  characterized  in  the  gas  phase.  Following 
Walsh,62  the  reaction  was  examined  for  a  linear 
configuration  along  the  reaction  coordinate  assuming  Cio 
symmetry,  with  the  energy  at  each  point  minimized  with 
respect  to  the  other  carbon-hydrogen  bond  distances  and 
angles.  The  barrier  for  reaction  (20)  is  0.52  eV  using  po¬ 
tential  I,  which  agrees  with  the  experimental  energy 
threshold  of  0.52  eV,63  and  which  can  be  compared  to  the 
value  of  0.58  eV  for  the  barrier  estimated  by  Walch  using 
ab  initio  total-energy  techniques.62  The  barrier  along  the 
linear  reaction  coordinate  for  potential  I  is  at  a  carbon- 
hydrogen  bond  distance  of  1.30  A  and  a  hydrogen- 
hydrogen  bond  distance  of  1.02  A.  The  values  calculated 
by  Walch^are  1.47  A  for  the  carbon-hydrogen  distance 
and  0.93  A  for  the  hydrogen-hydrogen  distance.  Poten¬ 
tial  II  yields  a  barrier  of  0.72  eV  at  a  carbon-hydrogen 
bond  distance  oJT  1.19  A  and  a  hydrogen-hydrogen  bond 
length  of  1.14  A.  This  barrier  is  somewhat  higher  than 
the  experimental  value,  and  so  potential  I  should  yield  a 
better  description  of  hydrogen  abstraction. 

IV.  E.  jERGETICS  FOR  STRUCTURES 
RELATED  TO  DIAMOND  CVD 

To  begin  to  quantify  energetics  for  surface  structures 
associated  with  the  CVD  of  diamond  films,  we  have  cal¬ 
culated  the  energetics  for  hydrogen  abstraction  from  the 
diamond  {111)  and  { 001  j  surfaces,  and  the  energetics  for 
the  chemisorption  of  methane,  acetylene,  and  their 
respective  radicals  at  terrace  sites  on  the  { 1 1 1 )  surface  of 
diamond.  More  systematic  studies  of  structures  and  dy¬ 
namics  of  hydrocarbons  reacting  with  diamond  surfaces 
and  their  relevance  to  the  CVD  of  diamond  films  are 
planned  to  be  presented  at  a  later  time. 

Removing  a  hydrogen  atom  from  a  hydrogen- 
terminated  fill)  surface  was  found  to  be  endoergic  by 
4.13  eV  for  potential  I  and  4.15  eV  for  potential  II. 
These  numbers  are  close  to  the  estimate  obtained  by  sub¬ 
tracting  the  relaxation  energy  given  in  Table  V  for  the  re¬ 
laxed  bulk-terminated  surfaces  from  the  carbon-hydrogen 
bond  strength  of  4.36  eV.  Combining  these  values  with 
the  gas-phase  bond  strength  of  4.75  eV  for  H2  yields  a  net 
exoergicity  of  =0.6  eV  for  abstraction  of  hydrogen  from 
the  {111)  surface  by  a  gas-phase  hydrogen  atom.  Be¬ 
cause  of  the  bonding  topology  of  this  surface,  the  energy 
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TABLE  VI.  Predicted  energetics  and  intramolecular  carbon-carbon  bond  lengths  for  various  single  molecules  chemisorbed  on  ter- 
race  sites  on  a  hydrogen-terminated  diamond  [111}  surface. 


Molecule 

Potential  I 
Potential  energy 
(eV) 

Bond  length 
(A) 

Potential  II 
Potential  energy 
(eV) 

Bond  length 
(A) 

Hydrogen  atom 

-4.1* 

-4.2' 

Methyl  radical 

-3.7* 

-4.0* 

Acetyl  radical 

-3.9* 

1.20 

-4.r 

1.29 

Hydrogen  molecule 

—  3.6b 

— 3.6b 

Acetylene 

—  5.0b 

1.33 

— 4.9b 

1.39 

Ethylene 

— 4.3b 

1.59 

— 4.3k 

1.57 

'Relative  to  a  hydrogen-terminated  surface  with  one  radical  site  and  the  gas-phase  molecule. 
bRelative  to  a  hydrogen-terminated  surface  with  two  adjacent  radical  sites  and  the  gas-phase  molecule. 


required  to  femove  a  second  adjacent  hydrogen  atom  is 
essentially  the  same  as  that  required  to  remove  the  first, 
so  radical  sites  on  this  surface  should  be  randomly  distri¬ 
buted.  Removing  a  hydrogen  atom  from  the  monohy¬ 
dride  phase  of  the  dimer  reconstructed  diamond  (001) 
surface  was  found  to  be  endoergic  by  4.44  and  4.49  eV  us¬ 
ing  potentials  I  and  II,  respectively,  so  abstraction  by  a 
gas-phase  hydrogen  atom  is  exoergic  by  ~0.3  eV.  Re¬ 
moving  a  second  hydrogen  atom  from  the  other  carbon  in 
the  surface  dimer  pair,  however,  is  endoergetic  by  only 
2.9  eV  with  potential  I  and  3.1  eV  with  potential  II.  This 
decrease  in  endoergicity  is  consistent  with  the  formation 
of  a  partial  double  bond  between  the  two  surface  carbon 
atoms.  Hence  the  dimer  reconstructed  (001)  surface 
with  radical  sites  should  be  composed  primarily  of  a  mix¬ 
ture  of  hydrogen-terminated  carbon-carbon  dimers  with 
single  bonds  and  nonhydrogen-terminated  double-bonded 
carbon-carbon  dimers. 

To  quantity  the  energetics  of  adsorption  on  terrace 
sites  on  the  (111)  surface,  chemisorption  energies  for 
both  single  molecules  and  monolayer  coverages  of  ace¬ 
tylene  and  ethylene  molecules  and  ethynyl  and  methyl 
radicals  have  been  calculated  using  the  empirical  poten¬ 
tials.  For  the  chemisorption  of  single  radicals,  a  hydro¬ 
gen  atom  on  a  hydrogen-terminated  (111)  surface  was 
replaced  with  the  molecule  (with  the  carbon-carbon  bond 
in  the  ethynyl  radical  perpendicular  to  the  surface),  and 
the  systems  were  relaxed  to  the  minimum-energy  struc¬ 


tures.  For  acetylene  and  ethylene,  two  hydrogen  atoms 
bonded  to  adjacent  surface  carbon  atoms  were  removed 
and  the  molecules  were  placed  with  the  carbon  atoms  in 
each  of  the  molecules  approximately  over  the  open  sites 
(i.e.,  with  the  intramolecular  carbon-carbon  bonds  paral¬ 
lel  to  the  surface).  The  resulting  energies  and  intramolec¬ 
ular  carbon-carbon  bond  lengths  for  each  of  the  mole¬ 
cules  are  given  in  Table  VI.  Energetics  for  hydrogen  ad¬ 
sorption  is  given  for  comparison.  The  radicals  CjH  and 
CHj  both  form  carbon-carbon  single  bonds  without  large 
changes  in  bonding  character  within  the  molecules, 
which.  le~ds  to  the  similar  binding  energies.  For  the  ad¬ 
sorption  of  acetylene  and  ethylene  each  molecule  forms 
two  single  bonds  with  the  surface  and  the  order  of  the 
carbon-carbon  bond  in  the  molecules  is  reduced  by  1. 
The  carbon-carbon  bond  in  the  chemisorbed  acetylene  is 
predicted  to  be  equal  to  the  gas-phase  double  bond 
length.  Hence  in  contrast  to  the  dimer  reconstruction  on 
the  (001 )  surface  the  constraint  of  the  lattice  does  not  in¬ 
hibit  the  formation  of  a  double  bond.  The  difference  in 
energy  between  the  chemisorption  of  acetylene  and 
ethylene  is  consistent  with  their  heats  of  hydrogenation. 
Hydrogenation  of  acetylene  to  ethylene  is  exoergic  by 
1.72  eV  relative  to  H2  plus  acetylene,  while  the  hydro¬ 
genation  of  ethylene  to  ethane  is  exoergic  by  1.35  eV. 

Given  in  Table  VII  are  energies  for  the  chemisorption 
of  a  monolayer  of  each  of  the  molecules  discussed  above 
on  the  (111)  surface.  In  each  case  the  initial 


TABLE  VII.  Predicted  energetics  and  intramolecular  carbon-carbon  bond  lengths  for  a  monolayer  of  various  molecules  chem- 
isorbed  on  the  diamond  (111]  surface.  The  energies  are  relative  to  a  relaxed  clean  surface  and  the  gas-phase  molecules. 


Molecule 

Potential  I 

Potential  energy 

Bond  length 
(A) 

Potential  11 

Potential  energy 

Bond  length 
(A) 

Hydrogen  atom 

—4.2  eV/atom* 

—4.3  eV/atom* 

Methyl  radical 

—  3.2  eV/molecule* 

—  3.0  eV/molecule* 

Acetyl  radical 

—4.2  eV/molecule' 

1.20 

—4.3  eV/molecule' 

1.29 

Hydrogen  molecule 

—  3.7  eV/mo!ecu!eb 

—  3.9  eV/moleculeb 

Acetylene 

—  5.2  eV/moleculeb 

1.33 

—  5.2  eV/moleculeb 

1.39 

Ethylene 

—  4.0  ev/mo!eculeb 

1.58 

—  3.9  eV/molecu!eb 

1.57 

‘One  surface  atom  per  chemisorbed  molecule. 
"Two  surface  atoms  per  chemisorbed  molecule. 
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configurations  were  generated  by  removing  the  hydrogen 
atoms  terminating  the  surface  bonds  and  translating  the 
chemisorbed  single  molecules  across  the  surface.  Each 
configuration  was  then  relaxed  to  the  minimum-energy 
structure.  The  energies  for  adsorbed  hydrogen,  ethynyl 
radicals,  and  acetylene  molecules  decrease  (i.e.,  become 
more  strongly  bound)  from  single  molecules  to  mono- 
layers,  while  because  of  the  interactions  between  non- 
bonded  hydrogen  atoms  the  energies  for  methyl  radicals 
and  ethylene  molecules  both  increase. 

While  a  specific  mode  of  CVD  growth  cannot  be  unam¬ 
biguously  deduced  from  these  energies,  they  do  suggest 
some  possible  trends.  First,  abstraction  of  a  single  hydro¬ 
gen  atom  from  a  ( 1 1 1  j  surface  should  be  favored  over 
abstraction  from  a  dimer  reconstructed  {001 }  surface. 
Abstraction  of  a  second  hydrogen  atom,  however,  should 
be  strongly  favored  on  the  reconstructed  (001  j  surface 
since  the  bonding  topology  of  the  surface  reconstruction 
allows  formation  of  double  bonds  on  the  surface.  Second, 
displacement  of  a  hydrogen  atom  from  a  hydrogen- 
terminated  diamond  { 1 1 1  j  surface  by  either  ethynyl  or 
methyl  radicals  is  not  strongly  energetically  favored. 
Ethylene  and  especially  acetylene,  however,  can  energeti¬ 
cally  displace  hydrogen  from  terrace  sites  provided  that  a 
gas-phase  hydrogen  molecule  is  formed  (rather  than  two 
H  atoms).  A  reaction  sequence  that  might  lead  to  this, 
however,  is  unclear,  although  a  sequence  similar  to  that 
proposed  by  Frenklach  and  Spear  may  be  possible.1  Fi¬ 
nally,  the  potential  suggests  that  adsorption  of  a  mono- 
layer  of  methyl  radicals  on  the  { 1 1 1 J  surface  is  not  an  en¬ 
ergetically  favorable  configuration  compared  to  hydrogen 
adsorption.  The  predicted  energies  given  here  for  a 
monolayer  of  methyl  are  uncertain  at  best  since  nonbond- 
ed  steric  interactions  between  hydrogen  atoms  have  not 
been  included  in  the  potential.  A  better  estimate  could 
probably  be  made  using  other  potentials  such  as 
Allinger’s  molecular  mechanics.32 

V.  SUMMARY 

As  a  first  step  toward  using  molecular-dynamics  simu¬ 
lation  techniques  to  study  the  CVD  of  diamond  films,  we 
have  developed  an  empirical  potential-energy  function 
that  can  model  intramolecular  chemical  bonding  in  a 
wide  range  of  hydrocarbon  molecules  as  well  as  diamond 
and  graphite  lattices.  The  analytic  function  is  a  highly 


parametrized  version  of  TersofPs  empirical-bond-order 
formalism  which  includes  terms  that  correct  for  an  in¬ 
herent  overbinding  of  radicals.  Nonlocal  effects  have 
been  incorporated  via  an  analytic  function  that  defines 
conjugation  based  on  the  coordination  of  carbon  atoms 
that  neighbor  carbon-carbon  bonds.  Because  only 
nearest-neighbor  interactions  are  incorporated,  the  func¬ 
tion  is  very  quickly  evaluated  and  can  therefore  be  used 
in  large-scale  molecular-dynamics  simulations.  Future 
refinements  to  the  potential  include  incorporating  bar¬ 
riers  for  rotation  around  carbon-carbon  bonds  and  non- 
bonded  interactions  such  as  van  der  Waals  forces. 

Atomization  energies  predicted  by  the  function  com¬ 
pare  well  to  experimentally  derived  energies  for  a  wide 
range  of  hydrocarbon  molecules.  The  function  predicts 
the  ir-bonded  chain  reconstruction  to  be  the  energetically 
favored  structure  on  a  clean  diamond  ( 1 1 1 )  surface  in 
agreement  with  local-density-functional  calculations,  and 
correctly  predicts  hydrogen  adsorption  on  a  bulk- 
terminated  surface  to  be  favored  over  the  ir-bonded  chain 
reconstruction.  The  structure  and  energetics  of  the  di¬ 
mer  reconstructed  diamond  I00))  surface  and  its 
monohydride  phase  were  also  ex:  ed  using  the  empiri¬ 

cal  potential. 

Energetics  for  a  limited  number  of  static  surface  struc¬ 
tures  that  are  relevant  to  the  CVD  of  diamond  films  were 
also  calculated  using  the  empirical  potential.  These 
structures  included  the  energies  required  to  abstract  hy¬ 
drogen-  atoms  from  diamond  (111)  and  dimer  recon¬ 
structed  jOOl)  surfaces,  and  energies  for  methyl  and 
ethynyl  radicals  and  acetylene  and  ethylene  molecules  ad¬ 
sorbed  on  terrace  sites  on  the  (111)  surface  of  diamond. 
More  detailed  and  systematic  studies  of  bonding  energies 
and  reactive  dynamics  using  this  potential  will  be 
presented  in  the  future. 
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